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We study a theoretical model of virtual scanning tunneling microscopy (VSTM),— a proposed 
application of interlayer tunneling in a bilayer system to locally probe two-dimensional electron 
systems (2DES) in semiconductor heterostructures. We focus on tunneling in the ballistic regime, 
where we show that the zero-bias anomaly is suppressed by extremely efficient screening. Since such 
an anomaly would complicate the interpretation of data from a VSTM, this result is encouraging 
for the ongoing experimental effort. 



I. INTRODUCTION 

The availability of increasingly clean low-density two- 
dimensional electronic systems (2DES) has allowed ac- 
cess to a regime in which electron-electron interactions 
play a major role.— >^ Evidence is accumulating from 
transport measurements that the physics of this regime is 
much richer than was previously appreciated (see Ref. 
and references therein). In particular, while much is un- 
derstood about the two limiting cases, = 1/ ^J na^ — > 
(Fermi liquid) and Vs oo (Wigner crystal), experi- 
ments on systems with intermediate values of rg = 10—30 
reveal a host of unanticipated anomaliesi^ (Here, n is 
the real density of doped electrons or holes and as is 
the effective Bohr radius, ub = fi.^e/m*e^, where m* is 
the effective mass, e is the electron charge, and e is the 
dielectric constant of the host semiconductor.) 

Experiments on 2DESs have been based mainly on 
transport measurements on large (micron to millime- 
ter scale) samples. Information on the local struc- 
ture of electronic states could powerfully elucidate the 
physics underlying these transport measurements, pos- 
sibly including the recently proposed "electronic mi- 
croemulsion phases."^ Momentum-space probes^ and 
finite-frequency probes^ have provided some important 
insights, but the residual spatial inhomogeneity in even 
the cleanest low-density 2DESs favors use of real-space 
probes. Over the past decade, important progress has 
been made in locally probing 2DESs, for example see 
Refs. iiHSliilliSlil- However, the challenge that a 
2DES is generally buried of order 100 nanometers deep 
in a heterostructure prevents use of powerful conven- 
tional techniques, such as scanning tunneling microscopy 
(STM), and a comparable position-sensitive scanning 
probe that can map out the local density of states at 
low energy is still lacking for such buried structures. An 
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ongoing effort to develop such a capability - termed vir- 
tual scanning tunneling microscopy or VSTM^ - is based 
on tunneling into a 2DES not from a scanned metal tip as 
in STM, but rather from a second "Probe" 2DES grown 
above the 2DES of interest (henceforth "Sample 2DES" ) , 
within the same heterostructure. Since the barrier be- 
tween the two 2DESs can be made very low by proper 
design of the layer structure, and since the Probe 2DES 
is not perfectly compressible, the barrier can be tuned 
at a particular location by applying a voltage to a sharp 
metal tip positioned above the heterostructure surface 
(Sec Fig[T|). Tunneling between the Probe and Sample 
2DESs is then strongly enhanced locally below the tip, 
and the location of enhanced tunneling can be scanned 
across the Sample 2DES by scanning the metal tip above 
the heterostructure. 

In this paper, we introduce a simple model of the 
VSTM - two parallel 2DESs connected by tunneling at a 
single point - and address the feasibility of the VSTM at 
its simplest level. For its intended purpose, the VSTM 
should ideally meet the following criteria: (i) there should 
be sufficient tunneling near zero bias to allow probing 
the low energy physics of interest, (ii) the tunneling rate 
should be sensitive to the local density of states at the 
location of tunneling. 

In the diffusive limit, the well known "zero-bias 
anomaly" occurs as a consequence of the inefficient 
screening of charge in 2D. Specifically, in 2D the con- 
ductivity has untis of velocity and hence the Coulomb 
energy E{t) ~ /R{t) associated with adding a charge 
to a 2D system decays with time in proportion to the 
screening radius R{t) ~ at, leading to an action that 
logarithmically diverges at small bias for tunneling into 
such a system. The result is a strong suppression of the 
tunneling rate near zero-biaSj ^^d^'-'^^d^ violating criterion 
(i). Moreover, the tunneling rate has a dominant contri- 
bution from long-distance physics, violating criterion (ii). 

The central result in the present paper is that in the 
clean limit, even in 2D, screening is sufficiently efficient 
that the tunneling action at zero-bias is finite and hence 
no zero-bias anomaly occurs. In this regime, the tun- 



ncling rate can be calculated pcrturbatively and is hence 
proportional to the local density of states. Our results 
indicate that using a VSTM to probe low energy local 
density of states should be feasible, if the 2DES of inter- 
est is clean enough. 

The outline of this paper is as follows: In section II, 
we present the model, which treats the tunneling elec- 
tron as a two state system and the remaining electrons 
in the 2DES that interact with the tunneling electron as 
the "bath" degree of freedom. In section III, we calcu- 
late the tunneling rate to lowest order in the tunneling 
matrix element. Finally, in section IV, we discuss the 
implications of our results. 



consists of two 2DESs characterized by 2D Fermi liquids 
with electron densities pi and p2, respectively, separated 
by a distance a. A voltage bias Vj^ias is applied across 
a single tunneling center at the origin (see Fig. [5]). We 
treat the tunneling electron as a two-state system repre- 
sented by = ±1 in the limit of a small bare tunneling 
matrix element A. The tunneling electron interacts with 
the density fluctuations of the 2DESs via a Coulomb in- 
teraction. 

In the ballistic transport limit, the action for this sys- 
tem is 

S[pi,p2,,<Jz] ^ Sa + Sa,p, (2.1) 




FIG. 1: A schematic depiction of the proposed VSTM setup. 
Figure courtesy of Adam Sciambi. 



II. THE MODEL 



Pi K a/2) uif^ 




FIG. 2: A tunneling electron in state = -f 1, upper layer 
(probe 2DES), interacting with charge density at position r in 
the same layer and charge density at position r* in the other 
layer (Sample 2DES). 

We consider a simple model that captures essential as- 
pects of the VSTM setup sketched in Fig. [T] Our model 



dhJ 



(2.2) 



(2^) 



(2^ 



ptKp + dtVp (2.3) 



where a = (1/2) < [1 + a^i^q, iw)], [1 - 0-^(9, >, p =< 
pi{q,iuj),p2{q,iuj) >, and 



K = 
V = 



Xr'(9,*^) V{q) 

Ui{q)^V{q) 

U2{q)-V{q) 



(2.4) 



Here Sa is the bare tunneling action in the absence of 
any interactions, while S'o-.p is the action for the rest of 
the ("bath") electrons, which we treat in the context of 
linear response theory: Xi denotes the density correla- 
tion function in each layer i, V{q) denotes inter-layer 
interaction and V is the coupling between the tunneling 
electron and the "bath" electrons through the intra-layer 
(U{q)) and inter-layer (1/(9)) Coulomb interaction. From 
here on, for simplicity, we restrict ourselves to the sym- 
metric case xr^('?i*'^) — X2^ {It^'^) = X~^('?i*'^) ^ind 
Ui {q) = U2 (q) = U{q), although the general case can be 
treated in an identical fashion. 

Since Eq. (|2.ip is quadratic, pi's can be readily inte- 
grated out to yield 



where 



So [(Tz 



doj 
2^ 



^VbiasCTz (iuj) - ^Actx (iw) -t- 5*0, (2.5) 



dLud^q[{Uiq)-V{q)Y 



(27r)3 



X~Hq,tuj) + Viq) 



(2.6) 



Here, the effects of correlations in the 2DESS are encoded 
in x(q, iuj). We treat the correlation effects through RPA 
in the rest of this paper; however, the form of the action 
in Eqs. (|2.5II2.6[) is more general. (The semiclassical re- 
sults of Levitov and Shyto^ti^ can be reproduced in this 
formalism if x is taken to be the susceptibility of a diffu- 
sive 2DES.) 
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III. PERTURBATIVE CALCULATION OF THE 
TUNNELING RATE 

A. Evaluation of the Action 

We evaluate the action Eqs. (|2.5ll2.6p using the RPA 
expression for x i^i the clean limit: 



(3.1) 



where the bare density correlator at zero temperature is 



Here, Vq = kp/irvp is the bare density of states per vol- 
ume at the Fermi surface. Furthermore, we assume that 
the distance a between the two 2DESs sets the short- 
est length scale and hence serves as the UV momentum 
cutoff, and hence we can make use of the approximate 

expression U {q)~V (q) ~ 27rae^ + O (^{\q\ a)^^ to lowest 

order in \q\a. This makes it possible to explicitly perform 
the q integral in Eq. (|2.6p to yield 



i/cj2 + {vpq)'^ 



(3.2) 



'5'off(T) = Sa 



'a^{iuj)\ k{lu), k{uj) 



1/7 



2tt 



Air^iyoVpil + a) 



'vpy 



2 2a log (^)(l + a)^ 



(1 



(3.3) 



where we define the dimensionless parameter a = 
2'Kae'^Vo and the UV frequency cutoff I/tq = vp/a. Ap- 
plying the analysis of Levitov and Shytov^ to this ef- 
fective action, it is easy to see that the accommodation 
time, as well as the action at zero-bias, are finite; there- 
fore, the tunnelling rate can be computed pcrturbativcly 
in A, as we do explicitly in the next section. 



B. Tunneling Rate 



ScS in Eq. (jlll A}i is of the same general form as for the 
"spin-boson" problem, in which the heat bath is treated 
as a collection of Harmonic oscillators. The heat bath of 
phonons is typically defined in terms of a spectral distri- 
bution function, J(ti^), which is simply the Hilbert trans- 
form of the kernel k{lo) in Eq. pil Ap : 



K (to) 



1 



duj' 



J {J). 



(3.4) 



Noting that the addition of a frequency independent con- 
stant to k(lS) results only in an (unimportant) additive 
correction to the ground-state energy, we find that 



where 

A = 



J(i^) = Au?e-^l^, (3.5) 
vp (1 + a) 



(3.6) 



The low frequency behavior J ~ is conventionally 
classifiedi^ as "supcrOhmic" for a; > 1 (the present case), 



where perturbation theory is applicable, "Ohmic" for 
X = 1 (which is obtained in the diffusive case), and "sub- 
Ohmic" for x < 1, in which cases non-perturbative meth- 
ods are necessary. 

With the spectral function in hand, we can straightfor- 
wardly calculate the tunneling rate following the steps of 
Ref. llSl . to second order in the tunneling matrix element 
to obtain 



(3.7) 



-h 



Here I\{x) is a modified Bessel function of the first kind. 
£ = (X £p has the dimension of energy and A and 
r2 are defined in Eq. (|3.6p . We note that the effect of 
Coulomb interaction enters the tunneling rate through 
the renormalized tunneling matrix element 



A = A exp 




where r. 



(l/n7r) 



1/2, 



(3.8) 



is the ratio of the Coulomb 



interaction energy to the kinetic energy, n is the elec- 
tron density, and as is the Bohr radius. Tunneling is 
suppressed for lower density, i.e. for larger Ts- The 
most notable feature of our results in Eq. p.7p and Eq. 
(|3.8p is the existence of an "elastic" term proportional to 
8 (Vbias/ft^^), which dominates the tunneling rate in the 
Vbias ^ limit. This term is absent in the Ohmic and 
subOhmic cases, due to the vanishing overlap (infrared 
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catastrophe) between the az = unperturbed ground 
states. This proves the existence of a finite tunneUng 
amphtude at zero bias. 

IV. SUMMARY 

In general, the zero-bias anomaly in tunneling into 
2DES's reflects the qualitative effects of the Coulomb 
interactions on the tunneling process. While these ef- 
fects are interesting in their own right, in the context 
of a VSTM they would represent a barrier to obtaining 
useful data. Through an explicit calculation, we have 
shown that in a ballistic system with efficient screening, 
the tunneling rate in the limit of zero bias is not en- 
tirely suppressed by Coulomb effects. This implies that 
the VSTM will indeed be capable of probing low energy 



physics of clean 2DESs through tunneling. We note that 
the main purpose of the current paper was a proof of 
principle, hence we limited ourselves to the simplest pos- 
sible application of the VSTM. There are many other 
systems to which the VSTM might be applied, where 
other considerations may be necessary, including tunnel- 
ing in a magnetic field and tunneling into a non-Fermi 
liquid. These issues will be subjects of future studies. 
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